Reentrant BCS-BEC crossover and a superfluid-insulator transition in optical lattices 
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We study thermodynamics of a two-species Feshbach-resonant atomic Fermi gas in a periodic po- 
tential, focusing in a deep optical potential where a tight binding model is applicable. We show that 
for more than half-filled band the gas exhibits a reentrant crossover with decreased detuning (in- 
creased attractive interaction), from a paired BCS superfluid to a Bose-Einstein condensate (EEC) 
of molecules of holes, back to the BCS superfluid, and finally to a conventional EEC of diatomic 
molecules. This behavior is associated with the non-monotonic dependence of the chemical potential 
on detuning and the concomitant Cooper-pair/molecular size, larger in the BCS and smaller in the 
EEC regimes. For a single filled band we find a quantum phase transition from a band insulator to 
a BCS-BEC superfluid, and map out the corresponding phase diagram. 



Feshbach resonances have become an essential experi- 
mental tool in the exploration of interacting degenerate 
atomic gases, allowing, for a realization of Fermi super- 
fluidity tunable from a weakly-paired BCS to a strongly- 
paired molecular regime [1-6J. 

Confinement of atomic gases in optical lattices is an- 
other powerful technique for realizing and tuning strong 
correlations [7 , allowing experimental investigation of a 
variety of lattice quantum many-body phenomena, such 
as the superfluid to Mott insulator transition [8 . 

Naturally, recent attention has focused on the rich 
combination of Feshbach-resonant gases in optical 
latticesjOHTT). Although considerable progress has re- 
cently been made ^l2nl8^ , a general theoretical descrip- 
tion of such system is challenging even at the two-body 
level as it involves a projection of the Feshbach resonance 
physics onto the eigenstates of the periodic potential. 
Optical lattice shifts bulk Feshbach resonances and in- 
duces new ones. 

Yet it is possible to argue [19l [20] that in a sufficiently 
deep optical lattice (bandwidth much smaller than a 
bandgap), even in the vicinity of a broad Feshbach reso- 
nance, the single band tight binding model with on-site 
attraction (attractive Hubbard model) /7(a), a function 
of a vacuum scattering length a is sufficient to describe 
the physics. Likewise, for a narrow Feshbach resonance 
in a deep lattice a tight binding model of open-channel 
fermions and closed-channel bosonic molecules coupled 
by an on-site interconversion term, Eq. ([T]) is appropri- 
ate. 

In this Letter, focusing on the deep lattice regime we 
establish that for a broad resonance, a two-component 
Fermi gas above half filling undergoes a BCS-BEC 
crossover with the molecular condensate in the BEC 
regime composed of holes; below half filling the crossover 
is to a conventional condensate of diatomic molecules. 
The BCS and BEC regimes are separated by an analog 
of a unitary point where the thermodynamic properties 
of the gas are given by a universal function of the band 
filling fraction. For a narrow resonance [4l[2T], while be- 
low half filling the behavior is qualitatively equivalent 



to its broad-resonance counterpart, above half filling the 
crossover is non-monotonic and reentrant. Upon decreas- 
ing the detuning the phenomenology crosses over from 
a BCS regime to a BEC regime of molecules composed 
of holes, back to a BCS regime and finally to a BEC 
regime of diatomic molecules. This rich behavior, illus- 
trated in Fig. ([l^) is associated with a corresponding non- 
monotonic dependence of the chemical potential on the 
Feshbach resonance detuning. Concomitantly, the size of 
a Cooper pair/molecule changes non-monotonically, with 
the BCS and BEC regimes respectively characterized by 
large and small pair sizes in units of the interatomic spac- 
ing. This resonant lattice phenomenology can in princi- 
ple be probed by measuring correlations in the gas after 
its free expansion [22] [23], through RF spectroscopy [24], 
compressibility, and via details of the atomic cloud's den- 
sity profile [7 , as illustrated in Fig. [5] 

Finally, we find that for a narrow resonance, a fully 
filled lattice (band filling of 2) exhibits a quantum phase 
transition between a band insulator and a paired super- 
fluid, absent for a broad resonance or in the absence of an 
optical lattice. This latter transition can be best detected 
in the "wedding cake" density profile, a layered structure 
of the gas when placed in an overall confining harmonic 
potential, with a shell of a band insulator sandwiched by 
an inner superfluid core and an outer superfluid shell. 

We now outline the derivation of these predictions. We 
study fermionic atoms in an optical lattice which contains 
N sites and M atoms (atom filling fraction n = M/N) 
for both the one- and two-channel models. The starting 
point is to consider the tight-binding two-channel Hamil- 
tonian 
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Here, cj^, Cr,cr are the creation and annihilation opera- 
tors of the fermionic atoms at lattice site r with spin a, 
6j, 6r are creation and annihilation operators of bosonic 
closed-channel molecules, ji is the chemical potential, g is 
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the coupling, is the "bare" detuning, t (t^) is the hop- 
ping matrix element for the atoms (molecules) and (r, r') 
denotes pairs of nearest neighbor sites. In the absence of 
interactions, the bosons and fermions are free particles, 
with the tight-binding dispersion of atoms given by 

ek = — 2t(cos kx + cos ky + cos kz)^ = — /i, (2) 

where the lattice constant is taken to be 1. Consider- 
able progress in understanding the thermodynamics can 
be obtained through a mean-field approximation where 
closed-channel molecular field is replaced by a classical 
field: 



6q ~ 5 ^q^O, 



(3) 



with 6q a Fourier transform of 6, 



Diagonalizing the re- 
sulting quadratic Hamiltonian and varying corresponding 
ground state energy with respect to and B* gives the 
number and gap equations: 
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with the integrals over the entire Brillouin zone (BZ). 
The solutions give /i and \B\ as a function of detuning 

We first analyze these equations in a limit of a broad 
resonance, corresponding to taking and g to infin- 
ity, while keeping their ratio U = /i'q finite. In this 
limit, the closed-channel molecules 6q can be adiabati- 
cally eliminated (integrated out) reducing the Hamilto- 
nian Eq. ([T]) to that of an attractive Hubbard model with 
the interaction strength U . 

In this limit the Eqs. Q and ([5| become 
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where A = g\B\ is finite, while \B\ goes to zero in the 
broad resonance limit. 

An important feature of these equations (by particle- 
hole symmetry that holds even beyond the mean-field 
approximation [2^) is that the sign of the chemical po- 
tential /i depends on whether the band filling fraction 
n is above or below 1 (corresponding to above or be- 
low half filling), with /i = for n = 1 independent of 
A and U. Mathematically, the latter is captured by 
/ d^k^kF{\^k\) = for any function F(|&|) for /i = 0. 
To see that the sign of /i is independent of A and there- 
fore [/, suppose n > 1. Then for A = (corresponding 



to U ^0), the right hand side of Eq. (|6| reduces to the 
Fermi-Dirac step-function, obviously giving /i > 0. Now, 
for a nonzero A, /i must remain positive since /a crossing 
zero at any A would imply n = 1 , contradicting the n > 1 
assumption. Similarly, if n < 1, then /i < 0, independent 
of U. 
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FIG. 1: The normalized chemical potential in, two-channel 
model (a) (for g/t — 20a/2) and one-channel model (b). Col- 
ored dashed lines correspond to different density n of the gas. 
The leftmost dashed line in (a) has n > 2. The dots on black 
lines correspond to the threshold of the molecule formation 
in the vacuum. 



Solving Eqs. (6|7) numerically gives the relation be- 
tween /i and U shown in Fig. ^p). We thus observe that 
in this broad resonance regime, for n < 1 the system 
undergoes a conventional BCS-BEC crossover, reaching 
the BEC regime for strong attractive interaction U (or 
equivalently for reduced detuning vq), where large nega- 
tive chemical potential gives (negative) half the binding 
energy of molecules formed by pairs of fermionic atoms. 
In contrast, forn > 1, as a refiection of particle- hole map- 
ping between n < 1 and n > 1 fillings, the ja > chemical 
potential actually grows with increased attractive inter- 
action, crossing above the top of the band for large at- 
tractive U. This corresponds to a BCS-BEC crossover to 
a BEC of molecules of two holes in the Fermi sea with 
the chemical potential tracking half of their binding en- 
ergy. At exactly half filling (n = 1) there is no BCS-BEC 
crossover, with chemical potential remaining pinned ex- 
actly at /i = for all [/, consistent with particle-hole 
symmetry of the Hubbard model at half-filling. 



At n = and n = 2 the Eqs. (6|7) can be straight- 
forwardly solved analytically. For example, at n = we 
expect that A = and /i < — 6t, dropping below the bot- 
tom of the band. Then Eq. (|6| is automatically satisfied 
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where C — Jg^ (^^J^s 3-cos(fc^)-cos(fe„)-cos(fe,) ~ 0.505. 
Writing the right hand side of the gap equation in this 
form allows us to expand the cosines around k = 
and extend the integral to infinity (valid as long as 
-At/(6t) - K 1) to find 
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This is the relation at the lowest density, n = 0, 

corresponding to half of the binding energy of the zero 
momentum molecule formed by two atoms. In order for 
the molecule to form in 3D, U must exceed a threshold 
attraction, > 4t/(7, as is known in the absence of lat- 
tice potential^ We note, however, that at a finite 
center of mass momentum at the Brillouin Zone bound- 
ary, a molecule can form at an arbitrarily weak attraction 
even in 3D [18^. 

Similarly by particle-hole symmetry, at n = 2 close to 
the top of the band /i(/7) is also given by Eq. (|9]). The 
full /i(/7, n) dependence for < n < 2 can be obtain 
through a straightforward numerical or an approximate 
analytical analysis of Eqs. (|6j [7|. However, it should 
be kept in mind that these are in themselves only valid 
within a mean field approximation, except at very small 
or very large U (at n 7^ 1). 

We observe that at a special value U = At/C the inter- 
action strength drops out of Eqs. ([6j [t]), corresponding 
to a divergent scattering length, just like at the unitary 
point of the BCS-BEC crossover in the absence of a lat- 
tice potential. Although these mean-field equations are 
approximate, this general feature holds more generally, 
with the chemical potential at a unitary point U = 
given by a universal function of particle density 
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with the property /(2 — n) = —f{n) dictated by the 
particle- hole symmetry. In the dilute n ^ limit, we 
recover the lattice-free result 



/(n) ^ -6 + ^ (37r^n 



2/3 



(11) 



where ^ is the Bertsch parameter [27]. ^ ~ 0.4[! 
The full function /(n) is not known but can be computed 
numerically. 

In contrast to the above broad resonance limit (one- 
channel model), where n is restricted to < n < 2, in 
the two-channel model, Eq. Q the filling fraction is no 
longer limited by two particles per site, as closed-channel 
bosonic states can accommodate an arbitrary number of 
fermionic atom pairs even if all states in the fermionic 



band are occupied. We next study the ground state be- 
havior of the two-channel model encoded in Eqs. (4|5), 
that in the narrow resonance limit, g^/n <^ Ep are quan- 
titatively trustworthy across the entire phase diagram 
and in the opposite broad resonance limit reduce to those 
of the one-channel model discussed above. 

Solving these equations numerically leads to the ijl — vq 
phase diagram illustrated in Fig. ([l^). Its many features 
can be understood analytically, particularly for a nar- 
row resonance. It displays a band insulator (BI) and a 
paired superfiuid (SF) phases, depending on the range of 
the chemical potential and detuning. For large positive 
detuning vq^ closed-channel molecules are separated by 
a large gap above the fermionic band, leading to a weak 
attractive interaction for the atoms. Thus, for a partially 
filled fermionic band, < n < 2, the chemical potential 
jiin) sits within the band — 6t < ji < Qt^ and the ground 
state is a weakly-paired BCS superfiuid. Increasing the 
filling to n = 2~ pushes the chemical potential to the 
top of the band, /i(2~) = 6t. Since the fermionic band 
is then full at n = 2, a further increase in n can only be 
accommodated by populating the closed-channel molec- 
ular state. For large vq^ the chemical potential therefore 
jumps from /J^{2~) = 6t to /i(2+) ^ which thus 

determines the lower and upper phase boundaries of the 
n = 2 BI. 

Reducing the detuning brings down the molecular 
state and leads to its hybridization with the pairs of the 
fermionic band states and a concomitant increases in the 
attractive interactions. Below half-filling, n < 1 this 
leads to a monotonic emptying of the fermionic band as 
the BCS superfiuid crosses over to the molecular BEC, 
familiar from a narrow resonance BCS-BEC crossover in 
the absence of the periodic potential [4 . 

The phase boundary between the paired superfiuid and 
n = BI (vacuum) can then be found exactly, as it cor- 
responds to a limit of a two-atom ground state, with 
/ic(^o) = m(^07^ = 0, 5 = 0). While for large positive 
i^o, the phase boundary /ic(^o) = — 6t follows the bottom 
of the band, for < di true stable molecular bound 
state (not just a resonance) peels off from the bottom 
of the band following half of the molecular binding en- 
ergy. To see this emerge from Eqs. ( |4|5[ ), we set n = 0, 
B = and note that for /i < —6t (below the bottom of 
the band) the number equation, Eq. Q is automatically 
satisfied. The gap equation then gives 



^0 -ttC + 2/i - ^ 



4t 



87rt3/2 



(12) 



2C/(4t) - m the depen- 
—{i^Q — ^0)^ and crosses 



For detuning just below = 
dence is quadratic /ic(^o) + 6t 

over to linear behavior /ic(^o) ^ ^(^0 — ^0)5 following 
the closed-channel level, as expected from the lattice-free 
analysis [41. 



Similar to broad resonance, Eq. (12) also depicts the 
phase boundary of n = 2. On this boundary the 
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threshold value of detuning below which the molecu- 
lar bound state of two holes can first appear is given 
by = g'^C/{4t) + 12t, with /ic = 6t, indicated by 
a dot in Fig. ([l^). Near and below this threshold 
point, ^/ jii — 6t dominates and leads to a lower branch 
that grows quadrat ically, corresponding to half the bind- 
ing energy of two holes. This reflects the particle-hole 
symmetry near this point and is consistent with ear- 
lier analysis of the broad resonance limit. The upper 
branch anticipated on general grounds, at large positive 
detuning asymptotes to the linearly growing solution, 
/^c(^o) ^ ^0/2, corresponding to the nearly decoupled 
closed-channel state that forms the upper phase bound- 
ary of the n = 2 BI. 

From further analysis of Eqs. (4j5), for n 7^ 2 we find 
the ground state is a superfluid for all with the chem- 
ical potential /i(n, uq) a smooth function of the filling n 
and detuning (i.e., exhibits a finite superfluid com- 
pressibility), as indicated by dashed curves in Fig. ([l^). 
For 1 < n < 2, /i(n, i^o) displays a non-monotonic de- 
pendence with uq^ that leads to a reentrant BCS-BEC 
crossover of holes at intermediate detuning and a con- 
ventional one of atoms at large negative detuning. 

In contrast, at n = 2 the system undergoes a quantum 
phase transition, shown on Fig. ([2|i), from a band insu- 
lator to a paired-hole superfluid as is lowered below a 

4 

critical value i^c = ^0 ~ bii-n'^t^ corresponding to the tip 

4 

of the band insulator lobe at jic = m(2, Vq) = 6t + 
in Fig. 

The onset of superfluid order close where B is small, 
can be studied analytical ly b y expanding the number 
and gap equations, Eqs. (4|5) in Taylor series in 
This gives that inside the SF phase below at n = 2, 
/i(2,z^o)/^ is a line of slope 3/4. Around the critical 
point, the superfluid order parameter displays the stan- 
dard mean-field onset, B ^ — ^o)^^^- The paired 
condensate (encoded in closed-channel molecular) den- 
sity \Bf as a function of detuning, for different atom 
densities n is illustrated in Fig. ([2|i), clearly revealing 
the BI-SF at n = 2. 

Because the fermionic n=2 BI is continuously con- 
nected to a bosonic Mott insulator at filling of one boson 
per site, beyond mean- field theory (valid for narrow reso- 
nance) we expect this particle-hole symmetric transition 
to be in the d-\-l dimensional xy university class. Among 
a variety of other probes, such as thermodynamics, den- 
sity profile, time of flight, and noise, the BI-SF transi- 
tion can be detected via compressibility, which enters the 
speed of sound for the superfluid mode, that we show van- 
ishes at the critical point and grows as (n — 2)^^^ away 
from it|20]. 

To make contact with cold-atoms experiments, the 
trap potential V{r) must be incorporated. This can be 
straightforwardly done through the local density approx- 
imation, /i /iiocai(^) = f^ — where the local chem- 
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ical potential /xiocai(^7 ^7 ^0) is maximum in the center 
of the cloud and drops off at its edges, and the global 
/i(A^, uq) is set by the total number of atoms TV in the 
gas. The resulting radial density profile p{r) is simply 
determined by a cut through the bulk phase diagram in 
Fig. ^) with p{r) = n{p\oca.\{r, N.uq)). For > Vc and 
average filling above 2 this predicts an n > 2 superfluid 
core, surrounded by a shell of an n = 2 band insulator 
(with requisite "wedding cake" plateau), that is further 
surrounded by a superfluid at n < 2. This is shown in 
Fig. ^p) for different detunings. 
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FIG. 2: (a) Paired condensate (closed-channel molecular) 
density \B\^ as a function of detuning at different total atom 
densities n. At total density n = 2, the molecular density 
vanishes above certain critical detuning, signaling a quantum 
phase transition into a band insulator, (b) Radial atomic den- 
sity p{r) of a trapped gas in the local density approximation. 
The plateau corresponds to an n = 2 band insulator shell, 
sandwiched by an inner (n > 2) and outer (n < 2) superflu- 
ids. 

To summarize, we studied an s-wave Feshbach resonant 
Fermi gas in a deep lattice potential faithfully modeled 
by a single band two-channel model. We showed that for 
above half lattice filling it exhibits an interesting reen- 
trant BCS-BEC crossover phenomenology of paired holes 
and atoms associated with the nonmonotonic dependence 
of the chemical potential on detuning. For a single filled 
band we find a quantum phase transition from an n = 2 
band insulator to a BCS-BEC superfluid, and map out 
the corresponding phase diagram. We expect that these 
predictions should be testable with current state of the 
art experiments on Feshbach-resonant Fermi gases in op- 
tical lattices. 
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